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We investigate the minimal performance, in terms of beam luminosity and detector size, of a 
neutrino factory to achieve a competitive physics reach for the determination of the mass hierarchy 
and the discovery of leptonic CP violation. We find that a low luminosity of 10 20 useful muon 
decays per year and 5 GeV muon energy aimed at a 10 kton magnetized liquid argon detector placed 
at 1300km from the source provides a good starting point. This result relies on Ova being large 
and assumes that the so-called platinum channel can be used effectively. We find that such a 
minimal facility would perform significantly better than phase I of the LBNE project and thus 
could constitute a reasonable step towards a full neutrino factory. 
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The recent discovery of #13 is a major step to- 
wards the completion of the leptonic mixing ma- 
trix. The first hints in global fits soon were followed 
by T2K [T], MINOS g] and Double-Chooz [3], and 
by now Daya Bay, RENO and Double-Chooz have 
observed this parameter at 7.7(7, 4.9cr and 2.9a, re- 
spectively [4-6 . The remaining major unknowns in 
the lepton sector are the Dirac CP violating phase, 
5 and the ordering of the neutrino mass eigenstates, 
sgn(Arn3 1 ). The large value of #13 admits many 
methods to determine the mass hierarchy (MH), 
some of which do not rely on long-baseline oscil- 
lation experiments, and therefore, we mostly will 
focus on the measurement of the CP phase. 

In order to measure CP violation (CPV), an ap- 
pearance experiment and consequently, a neutrino 
beam, is mandatory. The combination of T2K [7] 
and NOfA SB, may provide a hint for CPV, but the 
possibilities to obtain evidence at the 3cr level are 
rather limited [5] . Therefore a next generation neu- 
trino beam is needed. 

A multitude of experiments has been proposed in 
order to observe CPV in the leptonic sector, see for 
instance Ref. [10 . In the U.S. context, the long 
baseline neutrino experiment (LBNE) is the cur- 
rently preferred option. The present design con- 
sists of a 700 kW conventional neutrino beam aimed 
at a 10 kton liquid argon (LAr) detector placed at 
L = 1300 km from the source [TT]. Due to its long 
baseline and relatively large matter effects, LBNE 
would be able to determine the MH at 3er for a 75% 
of the parameter space. The CPV discovery po- 
tential is, however, limited due to a lack of statis- 
tics. An upgraded superbeam would obviously yield 
a much better physics potential [12]. However, su- 
perbeams are eventually limited by intrinsic back- 
grounds and systematic effects: the large flux uncer- 
tainties, combined with the inability to measure the 



final flavor cross sections at the near detector, in- 
troduce large systematical errors [13j [14] which are 
very difficult to control. 

For the determination of the CP phase a similar pre- 
cision to that achieved in the quark sector is only of- 
fered by a neutrino factory (NF) [T3J [TS] . In a NF a 
highly collimated beam of muon neutrinos and elec- 
tron antineutrinos is produced from muon decays 
in a storage ring with long straight sections |16| . 
The muons stem from pion decays and the pions, in 
turn, are produced by proton irradiation of a thick 
target. The present NF design parameters [T7] are 
10 21 useful muon decays per 10 7 seconds, aimed to 
a 100 kton magnetized iron detector (MIND) placed 
at 2000 km from the source, with a parent muon 
muon energy of 10 GeV. In order to form an in- 
tense muon beam for acceleration and storage, muon 
phase space cooling is required for this default con- 
figuration. In addition, the fact that the neutrinos 
in a NF are a tertiary beam implies significant pro- 
ton driver intensities; in this 4MW proton 
beam plus its associated target station. Thus, the 
challenges for an early implementation of a NF are 
clear: high proton beam power and muon cooling. 

These technical challenges are to be contrasted with 
the advantages of a NF - there are no intrinsic back- 
grounds and the absolute neutrino flux can be de- 
termined to much better than 1%. Furthermore, 
the presence of both muon and electron neutrinos 
in the beam does allow for a measurement of all 
final flavor cross sections at the near detector. A 
detailed analysis of the impact of systematic un- 
certainties in neutrino oscillation experiments has 
been recently performed in Ref. [2] , where the key 
systematics affecting the different types of facilities 
were identified. In particular, it was shown that the 
main sources of systematics affecting a NF are mat- 
ter density uncertainties (see also Ref. [T5] ) . On the 



2 




FIG. 1: CP violation (left panel) and mass hierarchy (right panel) discovery potentials. The statistical significance 
(for 1 d.o.f.) is shown as a function of the true value of <5, in the case when only the golden and disappearance 
channels are included in the analysis (blue lines, "gol+dis"), as well as when the platinum channel is also considered 
(red lines, "gol+dis+plat"). Dotted red lines show the results for the TASD detector, while solid lines refer to LAr. 
For reference, the results for the preferred reconfiguration option for LBNE are also shown (dashed gray lines). 



other hand, it was shown in Ref. [T5] that, for a 
large value of #13, the use of the so-called platinum 
channel, v ^ — > v e and — > v e , mitigates the effects 
from the matter density uncertainty. However, the 
platinum channel is inaccessible in a MIND, because 
it requires to identify electron neutrino charged cur- 
rent events. The feasibility of electron charge iden- 
tification was studied in Ref. [20] in the context of 
a low energy NF for a totally active scintillator de- 
tector (TASD) magnetized to 0.5 T using a so-called 
magnetic cavern. It soon was speculated that a mag- 
netized LAr detector should be suitable as well [21 
[23] . Note, that the detectors most likely will have 
to be deep underground due to the large duty factor 
of stored particle beams as in a NF. 

Specifically, we study a configuration with a muon 
energy of 5 GeV, 10 20 useful muon decays per year 
and polarity and a 10 kt magnetized LAr detector 
at a distance of 1 300 km. This choice of detector 
size and baseline is obviously inspired by LBNE: 
it allows to reuse the LBNE facilities to the largest 
possible extent and thus makes this option consider- 
ably more efficient in terms of resource usage. These 
parameters correspond to an overall reduction in ex- 
posure by a factor 100, and would in principle allow 
the NF to be based entirely on existing infrastruc- 
ture and technology. To give a specific example - 
the reduced beam luminosity would allow to use the 
existing Fermilab proton beam from the Main Injec- 
tor with a beam power of 700 kW and a beam energy 
of 60-120 GeV. This would result in a reduction of 



luminosity by a factor 10 with respect to the de- 
fault setup [17] - a factor ~ 5 from the beam power 
and factor 2 from the reduced pion production effi- 
ciency at higher beam energy. Furthermore, giving 
up muon cooling would result in the loss of another 
factor of 2 in luminosity. Conversely, an increase 
in the annual operation time from 10 7 s to 2 x 10 7 s, 
commonly used for Fermilab experiments, would re- 
sult in a gain of a factor of 2. As a result we finally 
obtain the postulated 10 20 useful muon decays per 
polarity and year. This is but one example of how 
to achieve this entry level luminosity. Obviously, 
future proton beam facilities may have much more 
favorable parameters resulting in corresponding in- 
creases in luminosity. 

Tab. [T] summarizes the detector parameters used in 
this work. A magnetized detector is required at a 
NF in order to disentangle the appearance and dis- 
appearance signals. In the absence of a detailed 
study of the performance of a magnetized LAr de- 
tector, we have followed Refs. [HI [55]. Since the 
LAr performance is indeed uncertain, we also eval- 
uate sensitivities using the performance parameters 
of a TASD, which have been studied in greater de- 
tail [20] , The same migration matrices and energy 
resolution as for the LAr detector have been consid- 
ered. Energy dependent efficiencies for the signal, 
following Ref. [22], have been used in this case, see 
Tab. [I] As for the backgrounds, we have considered 
the same rejection capabilities as for the LAr for 
muon detection, while for the electron detection we 
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have considered 10 times larger backgrounds. Fi- 
nally, 40 Ar cross sections have been used in both 
cases, which have been generated using GENIE [23] . 

In addition to the backgrounds considered in previ- 
ous references, the r-contamination [2STI27] has also 
been included in this work. In this case, no rejec- 
tion capability at all is considered when the polar- 
ity and flavor of the lepton produced in the r-decay 
is the same as that of the signal; while when the 
polarity is the opposite a 99.9% background rejec- 
tion has been assumed. Migration matrices, which 
account for the kinematics of the r-decay and the 
mis-reconstruction of the events at lower energies, 
have been generated using GENIE. A 17% r lep- 
tonic branching ratio has been assumed in all cases. 
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TABLE I: Main details used to simulate the LAr and 
TASD detector responses. The two rows correspond to 
the details used for and v e detection. The different 
columns indicate: signal efficiencies, background rejec- 
tion efficiencies (Neutral Current, Charge misID, Flavor 
misID), energy resolution and neutrino energy range. 
NC backgrounds have been migrated to lower energies 
using migration matrices produced by the LBNE collab- 
oration |12) . For details on the implementation of the 
r-contamination, see main text. 

Systematic uncertainties have been implemented as 
in Ref. [14], using the default values for the sys- 
tematic uncertainties listed in Tab. 2 therein. The 
only new systematic uncertainty introduced in this 
work is the one associated to v T cross sections, for 
which we have (conservatively) assumed a 45% un- 
certainty. All correlations are taken into account, 
between different channels as well as between the 
near and far detectors. All simulations have been 
done using a modification of GLoBES [28, 29], and 
marginalization over not-shown oscillation parame- 
ters was performed as explained in Ref. [14 . 

Figure [T] shows the results for the CPV and MH 
discovery potentials of the facility. These are de- 
fined as the ability of the experiment to rule out 
the CP conservation hypothesis and the wrong mass 
hierarchy, respectively, and their statistical signifi- 
cance is shown as a function of the true value of S. 
For reference, we also show the results for phase 



I of the LBNE experiment, which has been sim- 
ulated according to the CDR from October 2012, 
Ref. [TT]. Therefore, no r-contamination has been 
included in this case, and systematic uncertainties 
have been implemented as overall normalization er- 
rors over all signal and background contributions at 
once (no near detector has been simulated for this 
setup). Clearly, the low energy, low luminosity NF 
outperforms LBNE by a considerable margin once 
the platinum channel is included. 

The left panel in Fig. [2] shows the allowed region 
in 613-6 plane for one particular point in the pa- 
rameter space, where the different line styles corre- 
spond to different combinations of channels and/or 
experiments. Clearly, the addition of the platinum 
channel improves the performance beyond a mere 
increase of statistics - a true synergy, whose origin 
is explained in Ref. [19] . The right hand panel, on 
the other hand, shows the achievable precision for 
a measurement of S at 1 a as a function of the true 
value of 5. Again, we find that the low luminosity 
low energy NF constitutes a marked improvement 
over LBNE. 

We also would like to point out, that once a 4MW 
lower energy proton beam becomes available from 
Project X, muon cooling is added and the detector 
mass is increased by a factor 1-3, the performance of 
this facility matches or even exceeds the one of the 
baseline NF. Therefore, neither the initial energy of 
5 GeV nor the baseline do need to be changed in 
later stages. 

In summary, using 10 times less useful muon decays 
and a 10 times smaller detector mass with respect to 
the baseline NF still allows a muon decay based neu- 
trino beam to outperform realistic superbeam se- 
tups like LBNE, while offering at the same time the 
path to a full scale NF. Such low luminosity can be 
achieved using existing proton drivers and without 
muon cooling. Therefore, the timescale for an ex- 
periment like this would in principle be comparable 
to that of a superbeam experiment. Consequently, 
it lends itself well as an option for a second phase 
of an initial superbeam experiment. The reduced 
muon energy makes it possible to use a shorter base- 
line around 1300 km, and the use of a magnetized 
LAr detector allows to fully exploit the physics po- 
tential of the platinum channel, which is crucial for 
the overall performance of the facility. At the same 
time, these choice for baseline and detector technol- 
ogy ensure a good synergy with the first stage of a 
superbeam program. The big open question is how 
well can a LAr detector perform this task. 
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FIG. 2: Left panel: confidence region in the #13 — 5 plane for a particular point in the parameter space, at la CL 
(2 d.o.f.). Right panel: precision achievable for a measurement in S at the la CL (1 d.o.f.), as a function of the true 
value of S. Results are shown when only the golden and disappearance channels are included in the analysis (blue 
lines, "gol+dis"), as well as when the platinum channel is also considered (red lines, "gol+dis+plat" ) . Dotted red 
lines show the results for a TASD detector, while solid lines refer to LAr. For reference, the results for the preferred 
reconfiguration option for LBNE are also shown (dashed gray lines). 
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